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ABSTRACT 


Equations of motion are derived for"the transient response, to 
shock wave, of a submerged shell with internal structures. A substructuring 
procedure, which does not require calculation of a system stiffness matrix, 
employed to obtain these equations in a general manner for arbitrary 
internal structures approximated bv finite elements. 












TABLE OF CONTENTS 


I INTRODUCTION . , 

II GENERAL KINEMATICS . , 

III EQUATIONS OF MOTION OF SHELL . . . . 

IV EQUATIONS OF MOTION OF SUBSTRUCTURE 

V CONCLUDING REMARKS . 

REFERENCES . 

FIGURES . 

APPENDIX . 

















INTRODUCTION 


In an earlier paper (Ref. [1]), equations of motion for the transient 
response, to a shock wave, of a submerged shell with internally attached 
structures were obtained by a substructuring procedure. Several quantities 
in these equations, relating to the attachments, were only broadly defined, 
and determined in detail only in an ad hoc manner for specific appendages. 

In this report, matrix expressions for these quantities are defined which 
permit their computation for arbitrary internal structures approximated by 
finite elements. 

The present report addresses the same basic problem as Ref. [1], but 
does not consider the details of the fluid loading on the shell, which may be 
found in Ref. [1]. The effects of external loads applied to the internal 
structures, which were not considered in Ref. [1], are included here. 

The substructuring procedure described herein is intended for use in 
transient response problems. Unlike many other substructuring schemes, which 
are used primarily to solve steady-state problems (see, e.g.. Ref. [2]), the 
present method does not involve the modes and natural frequencies of the combined 
system (shell and appendages). In addition, a system stiffness matrix is not 
required, since the present technique employs the interaction forces and 
moments at the junctions of the shell and the attached structures. 






II GENERAL KINEMATICS 


Consider a shell S of arbitrary geometry (see Fig. 1). Let n, 2 and 3 
denote mutually perpendicular directions, respectively normal and tangent to 
the middle surface of the shell. If d^, d^ and d^ denote corresponding 
displacement components of this surface, a 3 by 1 displacement matrix d(s, t) 
may be defined, in which s denotes an arbitrary point on the middle surface, 
t is time, and underlining indicates a matrix; d^ is assumed positive inward. 
Corresponding to this displacement matrix, the 3 by 1 jth shell mode 


r. ^ T., 

and a generalized mass 


Wj = I “ij 


dA 


in which m denotes the mass per unit of surface area A of the middle surface 
of the shell and the superscript T denotes a transpose. The shell modes, 
which satisfy the orthogonality condition 

I m ^ dA = 0 if k j 
A 

may be used as expansion functions for the shell displacements, l.e.. 


d(s, t) = I qg^(t) i^(s) 
k=l 


in which the are generalized coordinates for the shell and is 

the number of modes of the shell considered. 

Assume that an elastic structure o, entirely within the shell S, is attached 
to S at points Sj(I = 1, 2,..., N^), as shown in Fig. 2. Let o be approximated 
by finite elements with Interface (attached) nodal points I, coinciding with 
Sj, and Interior (unattached) nodal points a(a “ 1, 2,...,Ny), corresponding 
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to material points of a not In contact with S. The finite element representation 
of O Is assumed to Include rotational degrees of freedom. 

At each Interface point I, certain components of displacement and/or 
rotation of O must match corresponding quantities of S, depending on the 
type of support. Thus, at an Interface point I, define the vector 



where the upper partition, an by 1 matrix, contains those physical degrees 
of freedom of a which must match corresponding physical degrees of freedom of 
S, and the lower partition, an by 1 matrix, contains the remaining physical 
degrees of freedom of a active at Interface point I.*^ 

Since the interior points of a are not connected to S, define the vector 
of physical degrees of freedom of a at an interior point a as 


U 

-aa 


X 


-oa 


(5) 


It may be observed in Eqs. (4) and (5), and in what follows, that the symbol - 
above a quantity Indicates a quantity associated with physical degrees of 
freedom of O constrained to move with S, while the symbol above a quantity 
indicates a quantity associated with physical degrees of freedom of a not 

^ At a simple support, e.g., the displacements of o and S must be equal, 
but the rotations ofo and S are not equal. Thus, assuming two active 
physical degrees of freedom, the displacement of o is located in the 
upper partition of Eq. (4), and the rotation of o is located in the 
lower partition of Eq. (4). 







constrained to move with S. 


If the appropriate definitions of Eqs. (4) and (5) are applied to all 
nodal points of O, the physical degrees of freedom of a may be arranged in 
the form I 

X, = 


in which the upper partition, an N by 1 matrix, contains the terms from 
all interface points [Eq. (4)] and the lower partition, an N by 1 matrix, 
contains the 3^^. terms from all interface points [Eq. (4)] and the 3^ 
terms from all interior points [Eq. (5)]. 

In the sequel, the elements in the upper partition of Eq. (6) will be 
referred to as "constrained physical degrees of freedom of o", while the 
elements in the lower partition will be referred to as "unconstrained physical 
degrees of freedom of a". The matrix 3^ has size N by 1, where 

N = N + N 

is the total number of active physical degrees of freedom of a. 

From Eq. (3), the shell displacements at an interface point 1 may be 
written as 

d(s,, t) E d . I , (t) 4(s,) E 

k=l 

where is a 3 by Vg matrix with elements ^(s^) and £g is a Vg by 1 matrix 
containing the generalized coordinates of the shell. Similarly, the rotations 
of the shell cross section at an Interface point 1 may be expressed as 


Si ^ ^ is 


(9) 







in which denotes the 3 by 1 matrix of rotations about the n, 2 and 3 

axes in the mode ^ and is a 3 by v_ matrix with elements (ti^(s ). 

*1C —1 5) —I 

Continuity of appropriate components of displacement and rotation of 
shell and substructure at an interface point I requires that 



( 10 ) 


where D^, an by 6 matrix, is a local rotational coordinate transformation 
which transforms quantities expressed in the S-coordinate system into the 
a-coordinate system. Each row of contains at most three non-zero terms. 
Use of Eqs. (8) and (9) permits Eq. (10) to be written in the form 


where is an by Vg matrix given by 



( 12 ) 


Application of Eq. (11) to all interface points, with reference to Eq. (6) 
for proper ordering, permits the continuity between shell and substructure 
to be written as 

4 = (13) 

in which C is an N by Vg matrix having the form 
~ __ 

The coordinate transformation array in Eq. (10) is used for convenience 
in presentation. Actually, only a 3 by 3 coordinate transformation array 
is required at each Interface point, the rows of this array applying to 
both displacements and rotations. 
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Equation (13) expresses the constrained physical degrees of freedom of a in 
terms of the free-free modes of S.*^ 

As suggested in Ref. fSl. the response of a may be taken as the super¬ 
position of its static response to the actual motion of its supports and 
the dynamic response with respect to fixed supports, i.e.**\ 

Xo-ioSo+sSo (15) 

In Eq. (15), is the N by fixed modal matrix, where denotes the 
number of fixed-base modes, and 3 ^ is the by 1 matrix of corresponding 
generalized coordinates. The N by N matrix is obtained by considering 
the static response of a to an arbitrary support motion (support "settlement") 
when the nodal loads corresponding to the unconstrained physical degrees of 
freedom vanish. In such a case, Eq. (15) becomes 

%>static = £ ^ 4 >static <16) 

*) " 

Equation (13) replaces Eq. (5) of Ref. [1]: 

^(Sj, t) = d(Sj, t) 

**) 

Equation (15) replaces Eq. (4) of Ref. [1]: 

^(s, t) = I ^ 1 ^®^ £(s, t) 







As shown in Ref. [2], the matrix g (sometimes termed the "constraint 
modes") may be obtained by considering the linear static "force-displacement" 


in which K is the N by N unconstrained stiffness matrix for the substructure 
and R is the corresponding vector of nodal forces and moments. For the 
determination of Eq. (17) may be partitioned to correspond to Eq. (6), 


-11 \ -12 


since the only nodal loads are the reactions. It then follows from the lower 
partition of Eq. (18) that 


, = (X ). 


where the superscript -1 indicates an inverse. A comparison of Eqs. (19) 
and (16) shows that 


Thus, the elements of g are the static displacements and rotations (including 
rigid-body motion), corresponding to the unconstrained physical degrees of 
freedom of o, which result from successive unit displacements and rotations, 
corresponding to the constrained physical degrees of freedom of a, at 
attachment points I when all other constrained physical degrees of freedom 
of O are set equal to zero. 


The use of Eq. (13) permits Eq. (15) to be written as 

Xj - io So + £ C gg 







As suggested In Ref. [2], Eqs. (13) and (21) may be combined Into the single 
matrix equation 

-iS 

in which X is given by Eq. (6), the transformation matrix has the form 


and the array of generalized coordinates is given by 


In Eq. (22), X is an N by 1 matrix, £ is a v by 1 matrix, with 


and T is an N by V matrix. Equation (22) may be written in subscript 
notation as 

Oj ■ '• 2. N) 










Ill EQUATIONS OF MOTION OF SHELL 

The kinetic and potential energies of the combination of shell S and 
internal structure o may be written as 


" ■ I “J’SJ * 2 Xcj 


, S 2 2 1 N N 


in which the are given by Eq. (26), the are the elements of the 
partitioned stiffness matrix of Eqs. (17) and (18), and the are the 
elements of the partitioned mass matrix of a from an assumed lumped mass 
*) 

formulation . The mass matrix must be partitioned to correspond to Eq. ( 6 ), 

l.e., 

M = 



(29) 


where the upper left partition, an N by N matrix, and the lower right 
partition, an N by N matrix, are diagonal. 


The equations of motion of the shell may be derived from Lagrange's 
equations: 


dt^3q„/ 


9V 


(j = 1,2,...,v„) 


(30) 


If Eq. (26) is substituted into Eqs. (27) and (28), one 
equations which replace Eqs. ( 6 ) and (7) of Ref. [1]: 




obtains th( 







in which the are the generalized forces corresponding to the total fluid 

pressure p acting on the surface of the sJiell (see Fig. 1) and the Q are 
T S j 

the generalized forces from the action of 0 on S, exclusive of the contribut 

derivable from the kinetic and potential energies. As shown in the Appendix 

the Q . are the elements of the matrix 
SJ 


Sc = C £ R 


where R, an N by 1 matrix, is the vector of known applied nodal loads acting 
on a (see Fig. 2). 


It follows from Eqs. (26) and (24) that 
3X„, 9X^, 


(J = 1, 2..... Vg; 

k = 1, 2,.. ., N) 


which may be seen to reference terms in the upper and lower left partitions 
of Eq. (23). Thus, the use of Eqs. (27), (28), and (32) leads to 


E'slr) ■ ''j ’s) * \ *0, 


3V 2 V C r 

= U-UiT q^. + y T, . y K, . X ^ 
9qs^ 3 3 ^Sj kj ki ai 


The substitution of Eqs. (31), (33), and (34) into Eq. (30), with reference 
to Eqs. (18), (23), and (29), results in 

"i ’s.1 V? Osi - Otj - '■) <j • '• 2 . V 


where the P are the elements of the vecti 












11 


whose addends are defined f 


P-c x„ + V fcl2 


p - C‘&‘ ISI22 x„ k ka + >122 


The bracketed term in Eq. (38) may be recognized as the physical 
equations of motion of a. Hence, the contribution of Eq. (38) vanishes, 
and the equations of motion of S become 


‘^Sj ‘^Tj % 


in which the generalized forces Q . , corresponding to the forces and moments 


exerted by o on S, are 


'=kJ<«k='akVi ■S.l V 


which result from the dynamic reactions 


'ok ■ -<«k 'ok * 'kl 'ol> 


exerted by 0 on S 


Using Eqs. (22) to (24), the equations of motion of S may now be written 


Equations (40) and (41) are generalizations of specific examples 
given by Eqs. (17), (18); (22), (23); (Al), (A2); and (A12), (A13). 
(A16) of Ref. [1]. 
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* “j"! ‘'sj ' 


•’tj - I, h)h <J ■ '• ^.V 


where 6., is the Kronecker delta and the F., the elastic forces and moments 
Jk 1 

exerted by o on S and vice-versa, are the elements of the N by 1 matrix 

I = (kii + ki2 ^ hi ^ 

The effect of inertia forces appears in the double summation in Eq. (42) 

It may be observed that the equations of motion of the shell have stiffness 
coupling only between shell and substructure. 

In the equations of motion of the shell, the contributions of the inertia 
forces and the elastic reactions of O on S may be given a simple physical 
interpretation by referring to Fig. 3. Figure 3a shows an Idealized finite 
element model of a substructure O attached to a shell S, and Figs, 3b and 3c 
show the manner in which Eqs. (42) and (43) treat the connection of a to S. 

The effect of inertia forces results from the rigidly attached point mass of 
Fig. 3b, while the effect of the elastic reactions of o on S results from the 
"spring force" of Fig. 3c. 
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IV EQUATIONS OF MOTION OF SUBSTRUCTURE 


The equations of motion of the substructure a may be obtained fr< 


Lagrange's equations: 


(3 = 1. 2 ..... V^) 


In which the are the generalized forces corresponding to the known 


are the elements of the matrix 




It follows from Eqs. (26) and (24) that 


(j = 1, 2 ,..., V„; 


j+"’s k = 1. 2, 


which may be seen to reference terms in the upper and lower right partitions 
of Eq. (23). Since the upper right partition of Eq. (23) is a null matrix, 
the use of Eqs. (27), (28), and (46) leads to 

d . 3T . _ V A M V 
dt^3q^j^ Vn Wii 

The substitution of Eqs. (45), (47), and (48) into Eq. (44) results in 
'^akj^“k+N ^ak+N \+N,i ^ai " ° ° 


which may be written in matrix form as 


^ fo ] H22I I jg_r io fiSzi 1^22^ 










If Eqs. (22) to (24) are substituted into Eq. (50), and the matrix product! 


are carried out, there results 

I22 i i22 So + 2‘i21 ^^22 S) £ 3s 

(51) 

+ S 222 S ^ 3s - 5 


In view of Eq. (20), which 'lefines the constraint modes, the third 
term of Eq. (51) vanishes. Moreover, as discussed in Ref. [2], the fixed- 
base modes of o are the solutions of the eigenvalue problem 

i22%>3 • “aj !!22%>J 

where is the natural frequency of the jth mode, and hence are orthogonal. 

Thus, Eq. (51) may be simplified to yield the equations of motion of a : 

= (53) 

where and are diagonal matrices whose main diagonal terms are, 
respectively, the generalized masses and generalized stiffnesses of o, and 

G = 1522 £ (54) 

is a by N matrix whose columns are the expansion coefficients of a modal 
series representation of the constraint modes. In subscript notation, Eq, (53) 
becomes ^ fj 

’oj * %i %i ■ Vj h - “jk <5 ■ 1- 2.“a' <”) 

It may be observed that the equations of motion of the substructure o have 
Inertia coupling only between shell and substructure. 

*) 

Equation (55) may be recognized as a generalization of Eq. (A8) of Ref. [1]. 
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V CONCLUUING REMARKS 


The equations of motion for the transient response of a shell S with 
internal structure 0 are given by Eqs. (42), (43), (53), and (54). It is 
Important to note that no total system (S + o) modal stiffness matrix is 
required in the present formulation, as is required when eigenvalues and 
eigenvectors for the entire system are sought (see, e.g,. Ref. [2]). It 
should also be noted that, in view of the stiffness coupling betAireen S and a 
in Eqs. (42) and (43), the modal mass matrix to be Inverted in the solution 
of a transient response problem has size Vg by Vg, determined solely by the 
number of modes of the shell employed In the solution. Formulations with 
inertia coupling between shell and substructure require the size of the modal 
mass matrix to be inverted to increase as the total number of modes V increases. 

The physical stiffness and mass matrices of a, K and M of Eqs. (18) and 
(29), respectively, may be obtained by modifying a finite element computer 
code such as SAPIV (Ref. [4]). The constraint modes may be determined by 
means of Eq. (20), or by utilizing a solution option available with SAPIV. Hence, 
it is possible to evaluate all of the matrices required in Eqs. (42) and (43) 
by using only those rows of K and only those main diagonal terms of M which 
correspond to the constrained physical degrees of freedom of o. Once the 
expansion coefficients of the constraint modes [Eq. (54)] and the modal 
coefficients for the interaction forces and moments [Eq. (43)] have been 
evaluated, the constraint modes and the matrices K and M may be eliminated from 
the final form of the equations of motion of S and a. This elimination 
and the absence of a system modal stiffness matrix will reduce the 
computer core (central memory) required to solve a given problem. 





16 


The present formulation will also reduce the amount of computer 
time (central processor time) required to set up the matrices in the equations 
of motion, as compared with that required by a method employing a system modal 
stiffness matrix. By way of example, if the formulation of Ref. [2] is used 
for a transient response problem, the matrices 


must be evaluated. The computer time required to evaluate Eq. (56) may be 
excessive, especially if the problem under consideration requires many modes 
V and/or many physical degrees of freedom of a, and if K is not well banded. 

It should be pointed out that the rearrangement of physical degrees of freedom 
of a to produce the partitioning of Eq. (18) usually results in a poorly banded 
matrix. For a lumped mass formulation, the computer time required to evaluate 
Eq. (57) is not significant compared with that needed to evaluate Eq. (56). 

Thus, a reduction in the amount of computer time required to set up the 
equations of motion is obtained by the elimination of a system moda] 
stiffness matrix. A reduction in solution time is also quite probable. 

The free-free modes of the shell S may be determined by means of a computer 
code such as B0S0R4 (Ref. [5]), while SAPIV may, of course, be used to obtain 
the fixed-base modes of the substructure a . If a consistent mass formulation 
is desired for a, the equations of Sections III and IV may be modified In a 
straightforward manner. However, in such a case, SAPIV, which is based on 
a lumped mass formulation, must be replaced by a computer code employing a 
consistent mass formulation. Although the equations derived in Sections Ill 
and IV deal with one piece of Internal structure, the present method is easily 
extended to Include any number of internal appendages. 
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APPENDIX - GENERALIZED FORCES DUE TO APPLIED LOADS ON SUBSTRUCTURE 

The virtual work due to the known applied nodal loads on substructure a 
may be written as 

5W = (65^)^ R 

in which 6)^ is the N by 1 matrix of virtual displdcements and R is the 
N by 1 matrix of applied loads. In Ea. (Al), both the virtual displacements 
and the applied loads correspond to the unconstiained physical degrees of 
freedom of a . From Eq. (21), it follows that 

so that Eq. (Al) may be expressed as 


6W = 6g^(2 R) + 6^(c'^ R) 


The virtual work of Eq. (Al) may be written in terms of generalized 
forces as 

6W = 5 


(A2) 


where is given by Eq. (24) and 


5 = 


is a V by 1 matrix partitioned to correspond to Eq. (24). The use of 
Eqs. (24) and (A5) in Eq. (A4) results in 

6W . 2s + i 
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Since the variations of the generalized coordinates are arbitrary, a 
comparison of Eqs. (A3) and (A6) yields 


Ss = 5 

(A7) 


(A8) 


The components of generalized force associated with the generalized coordinates 
of the shell S are given by Eq. (A7), and those associated with the 
generalized coordinates of the substructure O are given by Eq. (A8). 
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